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We analyze the dynamical stability of black hole solutions in self-gravitating nonlinear electro- 
dynamics with respect to arbitrary linear fluctuations of the metric and the electromagnetic field. 
In particular, we derive simple conditions on the electromagnetic Lagrangian which imply linear 
stability in the domain of outer communication. We show that these conditions hold for several of 
the regular black hole solutions found by Ayon-Beato and Garcia. 



I. INTRODUCTION 



In the last few years, interesting exact solutions |l], |[ 0, ^, |^, |[ Q have been obtained in the context of nonlinear 
electrodynamics (NED) coupled to gravity. These solutions are interesting because for a suitable choice of the 
electromagnetic Lagrangian they describe static black hole spacetimes where the black hole does not exhibit the 
usual singularity at its center, but instead all fields are everywhere regular and fall off in the asymptotic region. The 
existence of these solutions does not contradict the singularity theorems || since the corresponding spacetimes are not 
globally hyperbolic and their stress-energy tensor does not satisfy the strong energy condition. Examples of regular 
black hole spacetimes have been given by Bardeen |J. While the corresponding Einstein tensor satisfies the weak 
energy condition these spacetimes do not satisfy the vacuum Einstein equations. As shown recently by Ayon-Beato 
and Garcia Q the non- vanishing Einstein tensor in the Bardeen model can be associated with the stress-energy tensor 
of a nonlinear electromagnetic Lagrangian and the Bardeen black holes can be obtained as exact solutions of an 
appropriate NED coupled to gravity. 

An important issue in field theories is the study of stability of a given static or stationary field configuration. 
Unstable solutions have usually less physical significance than stable ones since they are likely to decay to a stable 
configuration. For example, it is known that static spherically symmetric regular particle-like solutions in Einstein- 
Yang-Mills theories are linearly unstable and numerical studies suggest that (depending on the details of the initial 
perturbations) a perturbation either causes the system to explode or to collapse to a Schwarzschild black hole (which 
is linearly stable) 

In the context of regular black hole solutions in NED, we study the behavior of linear (classical) fluctuations of 
the gravitational and electromagnetic fields in the domain of outer communication. Basically, there are two possible 
outcomes: 1) The fluctuations grow exponentially in time which means that the background configuration is (at least 
linearly) unstable. In this case a natural question to ask is to what configuration (if any) the solution will decay 
to. 2) The fluctuations are bounded in time which means that the background configuration is linearly stable. This 
suggests that the perturbed configuration eventually settles down to the background configuration after the possible 
emission of gravitational and electromagnetic radiation. 

In this article, we derive simple conditions on the electromagnetic Lagrangian that guarantee the corresponding 
black hole solutions to be linearly stable. In order to do so we cast the perturbation equations into a system of coupled 
wave equations with symmetric potential. Linear stability can then be established if this potential is positive definite. 
We show that the latter condition holds if the electromagnetic Lagrangian and its first and second order derivatives 
satisfy appropriate inequalities. We check these inequalities for several of the regular black holes that have been 
found in self-gravitating NED and so prove that these solutions are linearly stable. However, we also argue that this 
situation is not generic. That is, one can easily obtain an unstable regular black hole solution by slightly changing 
the electromagnetic Lagrangian belonging to a regular black hole solution. 

A brief review of the theory of NED in general relativity and of the regular black hole solutions is presented in 
section II. In section III we derive the perturbation equations and cast them into the form of a wave equation with 
potential. We derive sufficient conditions for this potential to be positive and show in section IV that the latter 
condition implies the linear stability of the black holes. Conclusions are drawn in section [V. A demonstration of the 
fact that the Bardeen black holes are linearly stable is given in an appendix at the end of this article. 
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II. REGULAR BLACK HOLES IN SELF-GRAVITATING NONLINEAR ELECTRODYNAMICS 

In this section, we first review a few facts from NED and then focus on spherically symmetric solutions. In particular, 
we give a brief description of the regular black hole solutions obtained in El 0, H]. 

A. Equations of motion 

The action for NED in general relativity is 



S{g,A] = ±-J V=9d 4 x(jR-£(F) 



where R denotes the Ricci scalar with respect to the spacetime metric g^, F^ v — d^Av — d^A^ is the electromagnetic 
field, and C is an arbitrary function in the invariant F = iF^F^. In Einstein-Maxwell theory, C{F) = F, but here 
we consider more general choices of electromagnetic Lagrangians.[^3) The field equations are 

Gt = 2 [C F F» a F va - 61C] , (1) 
V„ (£ F F a ») = 0, (2) 
V p * F afl = 0, (3) 

where Hp = dC/dF and G^, V and * denote the components of the Einstein tensor, the covariant derivative and the 
Hodge dual with respect to g^ v , respectively. [24 We have included the Bianchi identities (||) for completeness and 



later use. 

It is also convenient to express the field equations in terms of the field P^ = CpF^ v and the Hamiltonian Ti. = 
2£ F F — C which we consider to be a function of P = \P tlv P ilv @: 

Gt - 2 [H P P» a P va - 6t (2H P P - H)\ , (4) 
V M P a ^ = 0, (5) 
V M {*H P P a n = 0. (6) 

However, one should notice that the set of field equations (called the F framework) are equivalent to the set 

of equations (0M0) (called the P framework) only if the map F i— ► P is invertible. Since P = CpF, this requires that 

dP 

/ — = Cp(C F + 2CppF) , and 
ob 

dF 

o 7^ ~prF> = n P (H P + mppp) . (7) 

oP 

Finally, we mention that there is a duality between the F framework and the P framework: Using the fact that 

LpHp = 1, L = 2H P P - H, *P" a * P va = P^ a P va - 2<P„P, 

it is not difficult to see that the equations can be obtained from the equations (|l],||,||) by performing the 

following transformation [^) 

g^v ^ g^v = g^v 

F^ ^ F I1V = *P I1V , (8) 

C{F) ^ C(F) ee -H{P) = -H(-F). 

Notice however that the FP duality transformation (||) is only a symmetry of the system if C(F) = C(F), i.e. if 
C{—F) + 2FLp(—F) = C(F). The only Lagrangian that satisfies this identity is C(F) = aF + A/2, with a and 
A two constants. If a = 1 this is the Maxwell Lagrangian with a cosmological term. In general, the FP duality 
connects theories belonging to different Lagrangians. (For a more detailed analysis of duality invariant nonlinear 
electrodynamics, see Refs. 0,0.) 
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B. Spherically symmetric solutions 



We now focus on spherically symmetric and purely electric solutions in the P framework, equations (^,|^,^). From 
this, one can also obtain purely magnetic solution in the F framework (but with a different Lagrangian) by applying 
the inverse of the FP duality transformation (^). In the remainder of this section, we assume that the Hamiltonian 
Ti-(P) is negative and satisfies the weak field limit 

H(P) ~ P, as P^O, 

although this condition will not be a crucial assumption in the stability analysis below. 

A spherically symmetric spacetime can locally be written as M x S 2 where M is a two-dimensional pseudo- 
Riemannian manifold. The metric takes the form 

g = g ab dx a dx b + r 2 dVt 2 , 

where dil 2 = d-d 2 + sin 2 •& dip 2 is the standard metric on the two-sphere S* 2 and g — g a b dx a dx b is a metric of signature 
(— 1, 1), where the indices a, b label some coordinates on the manifold M. It is not difficult to see that 

Pab = -^Vab , (9) 



with the other components of P^ being zero, satisfies the NED equations (gLph. Here, fj a b is the volume element 
corresponding to g and Q is a constant that represents the electric charge. 
Next, Einstein's field equations (Q) reduce to 

Gab = 2g a bH, 
GAb = 0, 
Gab = — 2<7abA 

where capital indices refer to the angular variables •& and (p. In Schwarzschild coordinates, 

g = -S 2 Ndt 2 + ^ 7 (10) 

these equations yield S = S(t) (which can be set to 1 by an appropriate redefinition of t) and 

d t m = 0, d r m = -r 2 H{P), P = -~%, (11) 
where the mass function m is defined by N — 1 — 2m/ r. In particular, it follows from (O) that spacetime is static 



(as long as N > 0), i.e. Birkhoff's theorem 
Given the Hamiltonian 7i(P), equations ( 



holds. 

Ill) can be integrated to yield the mass function m(r): 



m(r) = M+ / r 2 Hdr, (12) 



where M is an integration constant representing the ADM mass. As a consequence of the weak field limit, we have 
m(r) — > M — Q 2 /2r as r — > oo and the metric is asymptotically flat with the same asymptotic as the Reissner- 



Nordstrom metric. If the Hamiltonian Ti-(P) is such that the integral in equation (12|) converges as r — * 0, we can 
choose the ADM mass such that m(r = 0) = 0. In particular if 7i — * —A/2 as r — > where A is a positive constant, one 
obtains m(r) — > Ar 3 /6, that is, the metric approaches the one of de Sitter space as r approaches zero. Hamiltonians 
with a similar type of behavior at r = and r = oo have been constructed in Ref. [Q, ^[|. As an explicit example, 
we consider the one that has been used in Ref. : 

1 - 3x 9 / - x 5/2 



* (p) = f <iw-sMtt^ • <13) 

where x = ^ — 2Q 2 P = Q 2 /r 2 and s — \Q\/2m. The corresponding mass function is obtained from (|l2|): 

n M 1 Q 2 r 3 , r 3 
m(r) = M 



2(r 2 + Q 2 ) 2 (r 2 + Q 2 ) 3 / 2 
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Thus, if M = m, the metric is regular at r = and so is the only non-trivial component of the electric field 
E = F tr = TipPtr- 

E = Qr 4 ( - 5Q 2 15 M \ 

* V(^ 2 + Q 2 ) 4 2 (r 2 + Q 2 )7/2 y ■ 

Furthermore, if the charge to mass ratio (2s) is small enough, the function N(r) intersects the r-axis at some values 
Tct r h with r c < r/i, and r c representing the radius of an event and Cauchy horizon, respectively. Properties of the 
Penrose diagrams belonging to these solutions have been discussed in Q. As expected, the corresponding spacetimes 
possess the same global structure as the Reissner- Nordstrom black holes except that there is no singularity at r = 0; 
r = is the origin of the spherical coordinates. 

We end this section after some comments. First, it is clear that there is a wide class of Hamiltonians that have 
the qualitative properties described above, since one could deform slightly the Hamiltonian (^) without affecting 
its asymptotic properties. Next, as observed in || [13), regular electrically charged black holes are not well-defined 
when interpreted in the F framework. What happens is that since F = —E 2 /2 is negative but vanishes as P goes to 
zero or infinity, it must exhibit a minimum at some P* . Therefore, the condition (0) is violated at at least one point 
and the electromagnetic Lagrangian C{F) suffers from branching, see Ref. Next, while the stress-energy tensor 
corresponding to (O) satisfies the weak energy condition (because H. < 0) it violates the strong energy condition 
since spacetime approaches the one of de Sitter near the horizon. This, together with the fact that the spacetimes 
possess Cauchy horizons, explains why there is no contradiction to the singularity theorems. Finally, magnetically 
charged solutions can also be obtained by using the FP duality transformation (||), as mentioned above. An example 
of a magnetically charged solution is provided by the Bardeen model which has been interpreted as a self-gravitating 
nonlinear magnetic monopole in Q. 



III. PERTURBATION EQUATIONS 



In this section we investigate the stability of the solutions described in section II with respect to non-spherically 
symmetric dynamical linear fluctuations outside the event horizon. In order to do so, we use a recently introduced 
gauge-invariant perturbation formalism 0, [l5| and generalize the equations obtained by Moncrief ]l6| to nonlinear 
electromagnetic fields. These equations (which we shall call the pulsation equations) are a set of wave-like coupled 
equations with symmetric potential acting on gauge-invariant and unconstrained quantities which therefore capture 
the purely physical degrees of freedom associated to the linear fluctuations. Our perturbation formalism, which is 
described in more details in for vacuum perturbations and in 113] for the Einstein-Maxwell case, has the additional 
advantage to yield the pulsation equations for any spherically symmetric coordinates on the background. 

Notice that an alternative approach is given by the Newman-Pcnrosc formalism (the background has a metric of 
type D). However, we have not succeeded in decoupling the equations along the lines of |Q. As we will see at the end 
of this section, one does not necessarily have the equivalence between the spectra belonging to the pulsation operators 
in the odd- and even-parity sectors when the electromagnetic Lagrangian is nonlinear. This might be an indication 
that the Newman-Penrose approach is difficult when nonlinear matter fields are coupled to the metric. 

The reader who is not interested in the derivation of the pulsation equations may skip the remaining of this section. 



A. Preliminaries 

Our task is to linearize the equations around a spherically symmetric background, where we have to keep in 

mind that P^ — T-Lp^{d^A v — d w A tl ). In the following, we assume that the function Tip is regular and positive for 
all r > rh, where rh denotes the radius of the event horizon. 

In terms of small perturbations of the metric and the vector potential, Sg^ and SA^, respectively, the linearized 
equations become 

SG^ = 2 {HppP^P^ 6P + Hp (2P { «SP v)a - P'P/Sgtrp) 

-g^KHpSP-Sg^pHpP-H)}, (14) 

= V UP ai , - 26g a[a P^ + i g^Sg ap P a J , (15) 
where in the following the function 



k = 1 + 2Hp 1 H PP P 



(16) 
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will play an important role. Notice that in view of the inequality (|7|) and our assumption on Hp to be regular and 
positive, k is regular and non-zero outside the event horizon if and only if the portion of the background solution that 
lies outside the event horizon is well-defined in the F framework. 
The variation of P^ u can be expressed as 

5P»u = \ n 2 P C FF (P afi 5F a0 - HpP° a P p Jg ap ) P M „ + H^SF^, (17) 
where in terms of H.(P), we have 

Cpp = k ^Tip^Tipp . 

General expressions for the perturbed Einstein tensor are found, for instance, in p| . 

The perturbed quantities Sg^ and SA^ are subject to the following gauge transformations: With respect to an 
infinitesimal coordinate transformation 5x^ i— > 5x p + X^ , where X p is a vector field, 

5g^ i-> Sg^v + £xg^u , SA^ ^ SA^ + £ X A P , (18) 

where £x denotes the Lie derivative with respect to X p . Furthermore, the potential A^ is subject to U(l) gauge 
transformations generated by a scalar x- 

5A, t » 5A, t + d pX ■ (19) 

Because of the gauge invariance of the theory, the perturbation equations are invariant with respect to the transfor- 
mations ( |l8| , |l9| ). Since we want to separate these transformations from the purely physical degrees of freedom, our 
strategy is to introduce gauge-invariant combinations of the perturbation amplitudes and to recast all equations in 
terms of these amplitudes. 

Since the background is spherically symmetric, it is convenient to expand all perturbations in spherical tensor 
harmonics. Perturbations belonging to different angular momentum number £ and m and with different parities 
decouple from each other. We discuss only the sectors with I > 1 since the spherically symmetric sector (I = 0) has 
already been discussed in the previous section. 



B. Odd-parity sector 



In each sector tm linear perturbations with odd parity are parametrized by two scalar fields k and v and a one-form 
h = h a dx a on the two-dimensional manifold M : 

Sg a b = 0, 5g A b = h b SA, SgAB = 2W(a$b) , 
SA^dx^ = vS B dx B , 

where in terms of the standard spherical harmonics Y = Y tm the transverse vector harmonics Sb are defined by 
Sb = {*dY)B = f] A B \I aY ■ Here and in the following, quantities with a hat refer to the standard metric g — dfl 2 on 
S 2 and quantities with a tilde refer to the metric g on M. It is not difficult to show that v is invariant with respect 
to both coordinate and E/(l) gauge transformations. In the gravitational sector, we can construct the gauge-invariant 
one-form (see [L5[ ) 

hQ™) =h-r 2 d (J^) . (20) 

Since the perturbation equations are gauge invariant we can choose a gauge in which k = and h^ mv ^ — h to simplify 
the calculations and then replace h by h^ mv ' = h — r 2 d(r~ 2 k) in the final result. 



1. Gravito- electromagnetic equations 

Since in the odd-parity sector all scalar perturbations vanish, the expressions ( p^Jl5| , p7| ) simplify considerably in 
this case. For example we have SP — 0, and so 5P M „ = Tip^^d^SAy — d„5A^). Explicitly, 

SPpu = Up 1 [dv A S B dx B - £(£ + \)v YdQ] v , (21) 
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where we have used d(S B dx B ) = AYdtt = -£{£ + l)YdCl. 

Introducing equation (^l]) into formula (|l5|), using the background expressions (||) and the fact that g al3 dg a p = 0, 
the linearized electromagnetic equations yield 



J f {Hp 1 dv) 



£(£ + !] 



Hp^v + QM ( -r =0, 



where * and & = *d* denote the Hodge dual and the co-differential, respectively, with respect to the metric g.p^] 
Next, from expressions (14) and (pl| ) we find 

SG Ab dx b = 4 {Q(Uv) - r 2 Ch) |§ , 5G ab = 0, SG AB = 0. 

Using this in connection with the corresponding expressions for the perturbed Einstein tensor in the odd parity sector 
(seeRef. 0), 



vd\- 



8G Ab dx h = |c? 

SGab = -$hV( A S B ) , 
we end up with the manifestly gauge-invariant equations 



+ (^+l)-2 + r 2 Gl)^§i 



<? [Hp'dv] + ^X^ Hp^v = -Q*d 



where here and in the following A 



1)(£ + 2). Equation ( p5| ) is void if £ = 1 since in this case, VmSb) = 



(22) 

(23) 

(24) 
(25) 



2. Pulsation equations (£> 2) 



We first consider the case £ > 2. Then, equation ( p5| ) allows to express h^ mv ^ in terms of a potential "J, 

vA 



where the factors 1/y/X and r turn out to be convenient. Introducing also $ = W ^TLp 1 v, we can integrate equation 
( p3| ) and eventually obtain the following pulsation equations [p6[ : 

V4XHp~Q 



-A* 



" A 

















£(£ + 1) 



$ = 0, 



-np /2 AHp 1/2 + ^n P 



$ = o, 



where A = — d)d — V a V a is the Laplace operator with respect to the metric g. These equations have the form of a 
coupled system of wave equations with symmetric potential. We can rewrite them into the form 



(D + S) 



0. 



(26) 



where 



D 



S = 







-n P /2 w a Hp 1 w a np /2 I ' 



A 



\ZAXH p < 



V4\T-C P Q 



t{i + i) + ^-Hj 



As we will show in the next section in more detail, stability follows from the fact that the spatial part of the operator 
D and the matrix S are manifestly positive. 
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3. Pulsation equations (1=1) 
In the sector I = 1, A vanishes and equation (|23|) can be integrated to yield 




= 4Q z/ + 6 J, 



where J is an integration constant. Using this into equation (|2J) gives 

« P <? [H^dv] + (1 + 4H P ^Pj v = -m P ^ . (27) 

The general solution to this equation is the sum of a particular solution v\ and the general homogeneous solution. 
The particular solution, which is proportional to J, can be chosen to be a function of r only. Therefore, v\ represents 
a stationary axial excitation of the background configuration. For Tip = 1 an explicit solution is given by v\ = J/r 
which corresponds to the Kerr-Newman solution to linear order in the rotation parameter J/M |L4|]. For a NED, we 
have not found an explicit expression for v\ , but we anticipate that v\ may give rise to slowly rotating generalizations 
of the solutions found in jj], |[ || . The homogeneous solutions represent electromagnetic radiation which propagates 
in a stable way since the potential on the left-hand side of equation ( p7|) is positive (see also the next section). 



C. Even- parity sector 



Gravitational perturbations are more complicated in the even-parity sector than in the odd-parity one. The metric 
is now parametrized by a symmetric 2 tensor, H abl a one-form Q a and two scalar fields K and G in the following way: 

Sg a b = H ab Y, 

SfJaB = <3aV B Y~, 

Sg A B = Kg AB Y + Gr 2 (i7 A \7 B - ~ 9AB^j Y, 

where G is not present for 1 = 1 since ^V^Vb — \ qab^J Y = in this case. For I > 2 one can introduce the 
following gauge invariants Jl8| : 

H$ v) = H ab 2V (aPb) , =K- 2 V ^L + h> {l + i )G , 

where p a = Q a — (r 2 V a G)/2. With respect to an infinitesimal coordinate transformation parametrized by the vector 
field X = i a Yd a + fV A Yd A we have 

Pa^Pa+U, G^G+ 2/. (28) 

Since G and p a can both be set to zero by an appropriate coordinate transformation, we can assume that these 
quantities vanish when deriving the pulsation equations below. In this gauge, H^ b nv ^ = H ab , K^ mv " > = K and Q a and 
G vanish. Since the perturbation equations are gauge invariant, it is sufficient to replace H ab by H^ b and K by 
j^( mv ) m ^ ne fi na l expressions and so obtain the equations in an arbitrary gauge. 
The electromagnetic potential is parametrized according to 

SA^dx^ = aY + fidY, 

with a = a b dx b a one-form and \x a scalar on M. The corresponding electromagnetic field is 

SFpu = d^S A u - dJA^ = [da Y + {dn - a) A dY]^ , 

from which it is clear that the amplitude a = a — dfi is ?7(l)-mvariant. Using the background expression (^) and 
£xF = dixF = d(r~ 2 HpQrj ab X a dx b ) we find that with respect to an infinitesimal change of coordinates, 

a + H P Qi(Zadx a ). 

In view of the transformations ( |28| ) we see that we can construct a coordinate invariant one-form which reduces to a 

if Pa =0. 
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Electromagnetic equations 



Having discussed the gauge issues, we now proceed to the derivation of the perturbation equations. First, a small 
calculation using ([l7n yields 



SPn 



with 



7T 

7 



[inY + jAdY]^ , 



Hp 1 „,.k-1Q 
— — *da H 



2k r 2 



(29) 
(30) 



where k is defined in ([16]) and H — H a a . In terms of the scalar it and the one- form 7 we can write the electromagnetic 
equations as 



id 



r 2 ir+Q(H-2K) 



<?7 



(31) 
(32) 



Equation (|32|) is the integrability condition for equation ( 31 ) and can be used in order to introduce a scalar potential 
<j) according to 



7 = *c 



(33) 



Equation @ can then be integrated to obtain the algebraic relation |27| 

r 2 7r = -t{l + +—(H— 2K). 

Using the equations (^9|) and (|30| ) one then obtains the following pulsation equations for the electromagnetic field 

£(£+1) 



Hp 1 ^ {H P d4>) + K- 



^(2kK-H) = 0. 



(34) 



If we turn off the gravitational amplitudes (i.e. if H = K = 0), we see that the sign of the potential depends on the 
function n defined in (|l^). In particular, the equations are stable as long as k > 0. In the next section we will show 
that this is also a necessary condition for stability, i.e. we will show that if k < in an interval, the equations admit 
exponentially growing modes for sufficiently large I. 



2. Gravitational equations 

Since the electromagnetic field couples back to the gravitational field, we have to consider the linearized Einstein 
equations ( |l4| ) as well. The harmonic decomposition of the linearized Einstein tensor can be found in [ p"4"|| . The result 
is quite complicated, but the structure of the equations becomes somewhat more transparent if one splits the tensor 
H a b into its trace and traceless part and then introduces the one-form 

C = H T ab F (\J a r)dx\ 

where = H a fj — \g a bH denotes the traceless part of H^. The relevant components of the linearized Einstein 

tensor may then be expressed in the following way: 

(5G)Z = S (v a Vb - \&ab£\ Y, 
~g ab SG ab - {T + G ab Hj b F )Y, 
SG Ab dx b = (U + Ubg)\vaY, 

{8G) T a b F ^ a dx b = (v - l -G ab Hl F dr + G ab (C a - r -\7 a K)dx b \ Y, 
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where 



and 



Ubg = 2„ \G dr + (i) ac G b c y*dr) 



TF 
b > 



s 

T 
U 
V 



2 ' 

-<?C - \g{C, dr) + AK + -g(dK, dr) - \k - ^^H, 



N 



($C)dr + (idC)*dr 



<HK+ -H 



H- 



dr 



^C)^+ 1 -dg(C,dr) + e -^±^C 
r r 2r l 

+ ]-AK dr - dgtdK, dr) + ( Ar - 

2 \ 2r 



N 

dK + —dH. 

2r 



(35) 
(36) 
(37) 

(38) 



Here, N — (V a r)(V a r) agrees with the function N defined in the previous section, see equation (|Io|). Note that 
equation (^5j) is void when £ — 1 since then VaV_bF = — 3(]abY. 

In our case G a b = ^gabH is proportional to g a b, so we have G ab H^ = and Ubg — 0. Next, it is not difficult to 
see from the right-hand side of equation ( |l4| ) that SGab is proportional to qab so that 5 = 0. Then, it follows from 
equation (|35| ) that the trace of H a b vanishes, H = 0. (When £ = 1, H can be set to zero by an appropriate coordinate 
transformation). The remaining components of the right-hand side of (Q) yield the expressions 



T 



U = iHpH^dcj), 



V 



HrdK, 



(39) 

(40) 
(41) 



where equation (|33j) has been used in the derivation for U. Following the vacuum case |14|, |15j we first look at the 
component g(U, *dr): Using equations (|37j) and (flO)), we find 



*d[C - rdK] = -4?d 



Tip(j> -irdr 



where we have used the identity g(a,*dr) = — *(a A dr) for any one-form a. This motivates the introduction of a 
one-form Z defined by 

Z = C - rdK + AH P (t) ^dr. 

Since dZ = we can introduce a potential £ with Z — d(. Next, we look at the combination 2V — Tdr. Using 
||) and (UglJ), we find 



equations 



2rg ^Z-4H P <l>Qdr,drj 



+ £(£ + l)Z 



(A + 3(1 -N) + 2r 2 U) dK+[\ + 4H P ^j Kdr = 



Q 2 



(42) 



By virtue of the background equations we have 



d r [r (A + 3(1 - N) + 2r 2 H)] = A + 6d r rn + 6r 2 H + 2r 3 Hpd r P 



= X + AHp-^t 



and since Z — d(, we can integrate equation (142) to obtain|25 



Q 



2rg(d(, dr) + £(£ + 1)( + (rX + 6m + 2r 3 H) K - 8,NH P j = 0. 



(43) 
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Finally, we compute the combination rT + g(U,dr) which yields, using the equations (t 



and 



$Z--g(Z,dr) 
r 



K + A 



<? [H P %dr) +2NH P 



e(£ + i)Q 



(44) 



This equation may be used in order to express the scalar K in terms of the other perturbation amplitudes. 

To summarize, we have the perturbation equations (|34|) for the electromagnetic perturbation potential <j> and 
equation ( f43| ) for the gravitational perturbation potential £ where H = and K can be eliminated using equation 
(^||). The next task is to bring these equations in form of a wave-like equation with symmetric potential which is 
more suitable for a discussion of stability. 



3. Gravito-electromagnetic equations 

The pulsation equations can be brought into the form of a wave-like equation with symmetric potential after the 
following transformations: 

* = VA-^r, § = v^f^---rr > ), (45) 
a + A \ r a + A / 

where a = 3rAr — 4r 2 H = 6m/ V + 2r 2 H. For the moment, we restrict ourselves to the case £ > 2. Note that a is 
positive outside the horizon of a regular black hole (see below). 

In terms of the new potentials defined in (Ha), a long calculation yields the final pulsation equations 



* \ I Vii - ViX ™ p Q W \ ( ^ 

'22 



A UJ + ( t )1$ )=0 ' (46) 



where 



1 



£(£+l)\-2N\ + arAr 



2NXb 



r 2 (a + X) I y ' J r 2 (a + A) 

f I i" " T " I I -! II I. > ( ^ ~ 

V22 = K 



^ + 4^ (A + 1 - N - 2r 2 H + Wn) + K?'*tot)f + 

r 2 r 4 (a + X) v ) p p r 4 (a + A) 2 



1 , , , 2Nb 

W = r (A + 1 - JV - 2r 2 H + 2A^k) + — - , 

a + A v ; (a + A) 2 

and where we have defined 

O 2 

b = \ + 4H P \. 



Here, we have also used the background expressions 

rAr = rd r N = l-N + 2r 2 H, 
rHp l ${H P dr) = -1 - N - 2r 2 H + 2Nk. 

The equations (^6|) generalize the pulsation equations obtained by Moncrief (l(| to NED coupled to gravity. (For 
H = P the equations above reduce to the one obtained in 119].) 



4- Pulsation equations (£ > 2) 

We now proceed to show under which conditions the potential in the equation (^6|) is positive semi-definite. 

We first show that the function a — 6m/r + 2r 2 H. is positive outside an event horizon: Let r = r/, denote the radius 
of the event horizon. Then, m(rh) — r^/2 and < rhd r N(rh) = 2r 2 H.(r — ru) + 1. Therefore, a(r = Th) > 2. Next, 
using the field equations d r m = —r 2 TL, we have d r (ra) — 4H.pQ 2 /r 2 > 0. This shows that a cannot become negative 
or zero for r > ru- In order to prove that the potential is positive, we first absorb a term r -1 Ar appearing in V\\ and 
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the term ft p 1/2 AHp /2 appearing in V22 into the differential operator. We may then rewrite the pulsation equations 



into the form 



(D + S) 



0, 



(47) 



where 



D 



S = 



-n\ 



1/2 



r 2 (a+A) 
V4\Hp Q 
r 3 (a+A) 



Cl 



2Nb 
a+A 

1 2jVfc 
' o+A 





_ \/4AW7 



1/2 



r 3 (a+A) 
.1(1+1) 1 4-Hf 



2Nb 
a+A 



r 4 (a+A) 



C2 



2jVb 
a+A 



with 



ci = A - 

C2 = Ci 

iu = Cl 



1 - A 
-4Nk, 
-2Nk. 



2r 2 H, 



Assume that Tt < and k > for all r > r^. It then follows that a ,6 and Ci are positive since N < 1. Therefore, the 
diagonal elements of S are positive. In order to show that S is positive definite, it remains to check the sign of its 
determinant. A small calculation yields 



r 4 (q + A) 2 
A 

This is manifestly positive if 



O 2 

det S = k£(£ + 1) [(a + A)c x + 2 A A] + SNkHp-^- [£(£ + 1) - 2Nk] 



< 2Nk < 



(48) 



(if k = S is only positive semi-definite.) In particular, this condition is satisfied in linear electrodynamics since then 
k = 1. Provided that the weak field limit holds and that k is regular and positive near the event horizon, we have 
Nk — > 1 as r — > 00 and JVk — > as r — > r^, so in those limits the condition ( [48| ) is satisfied. 

5. Pulsation equations (1=1) 

For £ = 1 it does not make sense to rescale £ by -\/A since A = in this case. Considering $ = (/a instead of the ^ 
defined in (|45|), we see that the equation for $ decouples when £ = 1. Repeating the arguments above, it is easy to 
check that the corresponding potential is positive if k > 0. The equation for ^ can be rewritten as 



r 



r 2 V a 



V^HpQ 



2Nb 



w - 



(49) 



On the other hand, the perturbation amplitudes are not necessarily gauge invariant when t = 1. Indeed, as shown in 
Jl4[ , the potentials can be seen to transform according to 



C h-> C + ra /, 



with respect to an infinitesimal coordinate transformation parametrized by the vector field X = — r 2 V a fYd a + 
JV A Yd A where / satisfies V a (r 2 V Q /) = (This makes sure that we remain in a gauge with H = 0.) Therefore, 
\Jr 1— > * + r/ while the electromagnetic potential $ is gauge invariant. Having solved for $, equation ( [49| ) can be used 
in order to obtain the gauge-invariant part of ^. 



IV. STABILITY ANALYSIS 



In the previous section, we have shown that linear fluctuations around a spherically symmetric and purely electric 
solution are governed by a wave-like equation with symmetric potential of the form 



P\7 a p- 2 \7 a P 



S) u = 0, 



(50) 
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where V denotes the covariant derivative with respect to the metric g and where P is a positive definite symmetric 
matrix and S is a symmetric matrix. The vector-valued function u on M is a gauge-invariant combination of the 
perturbed metric and electromagnetic fields. We will show below that linear stability follows if the potential S is 
positive definite. In the previous section, we have shown that this is the case if the following conditions hold outside 
the event horizon: 

H<0, H P >0, (51) 
< Nk < 3, (52) 



where k = 1 + 2Hp 1 Ti.ppP. The first two conditions ( |5l| ) imply that the weak energy condition is satisfied outside 
the horizon, while condition ( p32] ) implies that the function F(P) is well defined outside the horizon, see (]?]). We have 
checked that these conditions are verified for the Hamiltonian ( |l3|) used in ffl], for several values of the parameter 
s. However, it does not seem to be difficult to modify the function Tt(P) slightly and to violate, for example, the 
condition k > in a small interval outside the event horizon. As we will show below, the background solution is 
unstable if k < at some point outside the event horizon. 

One way of analyzing the behavior of the solutions to equation fl50| ) is by introducing Schwarzschild coordinates 
(t, r), see equation (|io|), in which case (|5(]) can be written as 



(d 2 ~Pd r _ M P^ 2 d r ,P + NS)u = 0, (53) 

where d Tt = Nd r denotes the derivative with respect to the tortoise coordinate. (Notice that P does not depend on 
i.) One then considers perturbations which are smooth and vanish at the boundary points r — (r* = — oo) and 
r = oo (r* = oo). The spatial part of the wave operator in (|5^) is symmetric with respect to the L 2 product 



u,v) = I u'v dr* 



for such perturbations. Since that operator is also real, self-adjoint extensions exist (in this case the operator is 
self-adjoint on some Sobolev space), and the linear stability can be analyzed along the lines of M, EQ| using the tools 



from spectral theory. As a first step one can look at solutions of (53) with time dependency of the form u ~ e lu)t , in 
which case 

uj 2 u = {-Pd r ,P~ 2 d r ,P + NS) u. (54) 

If the operator on the right-hand side of (|54|) is positive, which is the case if S is positive definite, it follows that 
lu is real and positive and there are no unstable modes (if ui = there might be modes with linear growth in t). 
On the other hand, if S is negative the operator is not necessarily positive. If one can find an eigenfunction with 
negative eigenvalue, u> is complex and one has a solution that grows exponentially in time. While it is difficult to 
solve the eigenvalue problem explicitly in general, there are simple analytical methods to show the existence of at 
least one unstable mode pfj| . In fact, if one can find a function u with (u, u) < oo and (u, Hu) < 0, this establishes 
the existence of at least one negative eigenvalue and the background solution is linearly unstable. 

However, the analysis of the eigenvalue problem (|54|) only gives information on the point spectrum (or on quasi- 
normal frequencies, but the corresponding modes u do not satisfy (u,u) < oo). Generalizing the results obtained in 
|n| on should be able to show the following stronger stability result: If the potential NS is strictly positive for all 
r > rh and if u(t = 0, .), dtu(t = 0, .) are smooth and have compact support, there is a bound C < oo such that 

\u(t,u)\ <C 

for all t > and — oo < < oo. (If NS is nonnegative but fails to be strictly positive outside the event horizon, the 
perturbations might exhibit linear growth in time.) 

We now use the above arguments to show that if k < 0, the solution is unstable: Suppose that k\ = ft(ri) < 
for some r\ > r^. Assuming that k is continuous, there is a 5 > small enough such that n < k%/2 on the interval 
(ri — S, n + S) and n < on the interval (n — 26, n + 26) . Consider the test function 



(¥,$) = (y/iX-Wp Qf(r),rf(r) 



for equation fl47|) where the function f(r) is smooth, vanishes identically outside the interval (ri — 26, r\ + 28) and is 
identically one for r £ (ri — 6,r\ + 8). It is clear that (u,u) < oo, and u'Su = k£(£ + l)/ 2 . Furthermore, a small 
calculation gives 

(u, Hu) < AS fl + i) K 2 - 5\n 1 \l{l + I), 
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where K is a constant which is independent on I but large enough, such that 



Qrdr, (h\ 



/2 -1 



np /2 d r , (n P 1/2 rf 



< K. 



It is clear that (u, Hu) < for large enough angular momentum number I. Therefore, the background configuration 
is linearly unstable if k is negative somewhere outside the event horizon. 



CONCLUSION 



We have derived the pulsation equations governing linear fluctuations on a spherically symmetric and electrically 
charged black hole solution in NED coupled to gravity. Under a mild technical assumption (namely that Hp vanishes 
nowhere outside the event horizon) the pulsation equations in the odd-parity sector cannot exhibit modes that grow 
exponentially in time outside the event horizon. However, in the even-parity sector, the stability properties depend 
strongly on the Hamiltonian H.(P), and one does not necessarily have the equivalence between the spectra in the odd- 
and even-parity sectors via a supersymmetric transformation like for linear perturbations in Einstein-Maxwell theory 
||l7f . As a consequence, the quasi- normal frequencies might be different in the odd- and even-parity sectors if the 
electromagnetic Lagrangian is not linear. 

In terms of the dimensionless variable x = \/—2Q 2 P = Q 2 /r 2 , sufficient conditions for linear stability are 

(55) 
(56) 
(57) 
(58) 

for all < x < Xh (i.e. r > r/j), where Xh = Q 2 / r h corresponds to the value of x at the event horizon. While we have 
shown that a solution is unstable if Ti. xx > for some x < Xh (the case H xx > cor respon ding to a marginal case 
where growth which is linear in time might appear) we stress that the conditions ( ^B| , ^6 58 ) are only sufficient. For 
example, there might exist stable black holes that satisfy H xx < but violate the condition (|58|). 

The linear stability of the purely magnetic solutions considered in j|, ||] can be discussed on the same lines as above. 
In fact, using the FP duality transformation (||), one can easily translate the conditions ( |55| , p6| , p7| , p8| ) to obtain the 
sufficient conditions 



n 


< 


0, 


n x 


< 


0, 




< 


o, 


3H X 


< 


xNH 



C 


> 


0, 


r 


> 


0, 


r 

'-yy 


> 


o, 


3Cy 


> 


yNC 



yy ' 



(59) 
(60) 
(61) 
(62) 



where y ~ y /, 2g 2 F and g is the magnetic charge. 

One can show that the conditions ( |59|j60| , |6l| , |62| ) are satisfied for the Bardeen black holes of jl| for any value of the 
charge to mass ratio (such that one has a black hole). A proof of this statement is given in the appendix. We have 

I) hold for the black hole solutions found in 



also checked that the conditions (55 



p], ||, H for different values of 



the charge to mass ratio. This establishes the linear stability of these solutions with respect to linear fluctuations in 
the domain of outer communication. An interesting question is what happens to perturbations inside the black hole. 
As expected from the mass inflation scenario |2l[] fluctuations that are falling into the black hole are likely to destroy 
the Cauchy horizon and to introduce a singularity inside the hole. These questions will be analyzed further in some 
future work where we also plan to discuss the stability of particle-like solutions p2] . 
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APPENDIX A: LINEAR STABILITY OF THE BARDEEN BLACK HOLES 

The Bardeen black holes are obtained in the F framework from the Lagrangian Q 



n / \ 5/2 

3 / V • 



2sg 2 \l + y 



where y = y^2g 2 F and s — \g\/2m is half the monopole charge to mass ratio. Since 

15 y z ' 2 

L y 4 S£ ,2 (l+y)7/2 ' 



f(y) = yC-'C 



mi 



15 y 1 ' 2 
8sg 2 (1+y) 9 / 2 ' 
3-4y 
2(1 + 1/) 



it is clear that C y > and so it is sufficient to check the condition 

3 > Nf(y) > (Al) 

for < y < yh where yh is the value of y at the event horizon. 

Since f(y) is monotonously decreasing with /(0) — 3/2 the first part of the inequality (|Al|) is verified. Next, 
consider the component 

1 yl/2 

N = - 9tt = 1 -:s(i+ y f/ 2 

of the metric which determines the location of t he event horizon. N has a single minimum at y = 1/2, so yh < 1/2. 



Therefore, the second part of the inequality (Al) is also verified for all < y < yh and we conclude that the Bardeen 
black holes are stable. 
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